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Abstract 

Wc find a continuum of extinction rates for solutions u(y,r) > of the fast diffusion 
equation u T = Au m in a subrange of exponents m € (0, 1). The equation is posed in 
R ra for times up to the extinction time T > 0. The rates take the form ||u( , ,t)|| 00 ~ 
(T — r) e for a whole interval of 8 > 0. These extinction rates depend explicitly on 
the spatial decay rates of initial data. 



1 Introduction 

We consider the Cauchy problem for the fast diffusion equation: 

( u T = A(u m /m), y G M n , r G (0, T), 



(1.1) 

u(y,0) = u (y) > 0, t/£ 



p II 



where m G (0, 1) and T > 0. The factor 1/m is not essential; it is inserted into the equation 
for normalization so that it can also be written as u T = V-(u m - 1 Vu). In that way, it is 
readily seen that the diffusion coefficient c(u) = u m ~ l — > oo as u — > if m < 1, hence the 
name Fast Diffusion Equation (but notice that c(u) — > as u — > oo). Furthermore, it is 
known that for m below a critical exponent m c = (n — 2)/n all solutions with initial data 
in some convenient space, like L p (M n ) with p = n(l — m)/2, extinguish in finite time. We 
will always work in this range, m < m c , and consider solutions which vanish in a finite 
time. The purpose of this paper is to study the rates of extinction of such solutions. Our 
main contribution is to provide a continuum of rates of extinction for fixed m. Technical 
reasons imply that m must be in the range < m < = (n — 4)/(n — 2), n > 5, for the 
construction to work. This restriction may be essential. 
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Let us review the state of the question from a broader perspective. The description of 
the asymptotic behaviour of the global in time solutions of (jl.ip as r — > oo for m > m c is 
a very active subject, and the study has been extended in recent times to the behaviour 
near extinction for m < m c , both in bounded domains or in the whole space. In the 
former case, the rate of decay for bounded solutions is universal, of the form \\u(-, t)\\oo = 
0((T — t) 1 /^ 1- " 1 )) when m > m s = (n — 2)/(n + 2), cf. [HE], but the question is more 
complicated when m < m s . 

In the case of the whole space, which is the one of interest here, the book [TU] contains a 
general description of the phenomenon of extinction, where it is explained that not only 
the occurrence of extinction depends on the size of the initial data, but also that different 
initial data may give rise to different extinction rates, even for the same extinction time; 
this may happen for all < m < m c . It is also proved, cf. [10] and quoted references, 
that the size of the initial data at infinity (the tail of uq) is very important in determining 
both the extinction time and the decay rates. 

Special attention has been given recently to particular classes of data that produce definite 
estimates. This happens in the case of data with the maximal decay rate compatible with 
extinction in finite time, which is 



u (y) ~ A\y\ 



2/(1 - m) 



;i.2) 



as \y\ — > oo. Note that [i < n for m < m c so these data are not integrable. Thus, the 
papers [5J [2j 0] are concerned with the stabilization as r — > T of general solutions towards 
some special self-similar solutions Ud,t known as the generalized Barenblatt solutions, 
given by the formula 



U D ,T{y,r) 



1 



R(r) n 



D + 



0(1 



m) 



y 



R(r) 



(1.3) 



where for m < m c we put R(t) := (T — r) 



0:-- 



1 



and 
1 



/' 



n(l — m) — 2 n(m c — m) 2{n — fi) 

Here T > (extinction time) and D > are free parameters. Note that R depends on T. 
It has been proved that the corresponding Barenblatt solutions with exponent m > m c 
play the role of the Gaussian solution of the linear diffusion equation in describing the 
asymptotic behaviour of a very wide class of nonnegative solutions, i.e., those with initial 
data in L 1 (M n ), cf. [11]. To some extent, the solutions (j 1 .3f) play a similar role for m < m c 
but their basin of attraction may be much smaller. This is precisely described in [5], with 
results on the basin of attraction of the family of generalized Barenblatt solutions; it 
establishes the optimal rates of convergence of the solutions of (jl.lj) towards a unique 
attracting limit state in that family. All of these solutions will have a decay rate near 
extinction of the form 'r)|| o = 0((T — t) h ^), and it is clear that n/3 > 1/(1 — m). 

A very interesting limit case occurs if we take D = in formula (II. 3D . and we find the 
singular solution 

U , T (y,r) :=K{T-TY' 2 \y\-^ (1.4) 
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whose attracting properties have not been studied. Note the value fe* = (2(n — fj,))^/ 2 . 

The question that we address here is the following: Can we obtain different decay rates 
near extinction for bounded data uo(y) that behave at infinity in first approximation like 
the singular solution, i.e., uo(y) ~ ? We will show that the answer is yes, and 

actually we will obtain a whole continuum of rates. 

Theorem 1.1 Let u > be a solution of Problem assume that 

(n — 4 \ 
, ) , (1-5) 
n — 2/ 

and let the initial function uq be continuous, bounded, and satisfy the conditions: 

0<u (y) < A\y\-» for all y^^ 

and 

Mv\~^ -ci\y\~ l < My) < A \v\~^ -aW 1 f° r \v\ > 1 

for some A, c±,C2 > 0, and 



1 1 



+ 2<l<L = ix+yj2(n-iJL). (1.6) 



Then the solution has complete extinction precisely at the time T = (A/k*) 1 m > 0, and 
there are positive constants K\ , if 2 such that for < r < T we have 

K\(T - r) e < |K, r)|U < K 2 {T - r) 6 , (1.7) 
where 9 = > 0, 7 = fazgKzizO , 

2{n—fi) ' 1 I— (i 



It is easy to check that under the above assumptions 9 covers an interval [9 m i n , 9 max ) with 
< 

9min ^ 9 max — fj,n/2,(n /x) — n/3. This is the precise range of extinction rates of these 
solutions, to be compared with the standard extinction rate (T — r) n ^ of the Barenblatt 
examples. 

As a precedent to this result, the existence of different rates was established in Theo- 
rem 7.4 of [10] for all m < m c by means of the construction of self-similar solutions of 
the form u(y, r) = (T — r) a f(y (T — t)^). In this way a whole interval (a, 00) is covered, 
which extends the scope of our present theorem. However, a (the anomalous exponent) is 
not explicit, we obtain only one solution for each time-decay rate and the dependence of 
a on the spatial behavior of the data is not analyzed. Theorem [TTT] clarifies these aspects, 
explaining the delicate relationship between both limits, \y\ — > 00 for uq and r — > T for 
u(y,r). 

The proof of the theorem needs techniques that are only natural after rescaling the 
problem. In fact, the rescaled problem allows us to formulate and prove a more precise 
result about the dependence of the rate on the tail of the data and the convergence 
of the spatial shapes. We devote the next section to the presentation of the rescaling 
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transformation, the resulting rescaled equation and the asymptotic convergence plus grow- 
up result in that context. Sections 3-5 will be concerned with proving the result for the 
rescaled problem. The last section is devoted to comments and open problems. 

Notations. Throughout the rest of the paper and unless mention to the contrary, we 
keep the conditions n > 5 and m < m*. The exponent m* also plays a big role in the 
asymptotic results of [SlElSj- We also keep the above symbols and variables. In particular, 
fj, = 2/(1 — m) so that m < m c means \i < n and m < m* means fi + 2 < n. 



2 The rescaled flow 

As we have just said, it is very convenient to rescale the flow and rewrite in self-similar 
variables by introducing the time-dependent change of variables 



1-m, f R(r)\ , B(l-m) y . , 

t:= —^{w>) and *~y EL -r- L wy (2 ' 1) 

with R as above, and the rescaled function 

v(x,t) := R(T) n u(y,T). (2.2) 

In these new variables, the generalized Barenblatt functions Uo,T{y^T) are transformed 
into generalized Barenblatt profiles Vd(x), which are stationary: 

V D (x) := (D+ \x\ 2 )^, x£R n . (2.3) 

If u is a solution to then v solves the rescaled fast diffusion equation 

v t = A(v m /m) +/iV • (id), t>0, i£l n , (2.4) 

which is a nonlinear Fokker-Planck equation (NLFP). We put as initial condition vq{x) := 
R(0)~ n uo(y), where x and y are related according to (|2.ip with r = 0, x = cy. Roughly 
speaking, vq is a rescaling of uq depending only on T. We have taken the precise form 
of this transformation from [3]. Note also that the factors 1/m and \jl in equation (|2.4p 
can be eliminated by manipulating the change of variables, but then the expression of 
the Barenblatt solutions would contain new constants. Thus, in our scaling the singular 
solution becomes 

V (x) = \x\~ IM , xeR n \{0}. (2.5) 



2.1 Main result for the NLFP equation 

In the following sections we consider the f-equation (|2.4p with initial data given by a 
bounded function < Vq < Vo, and such that the difference Vq — vq has a tail controlled 
by a power rate. This is our detailed result about asymptotic behaviour of the solution 
whose initial data vq(x) are perturbations of the steady state Vq(x). 
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Theorem 2.1 Assume that n and m are as in (jl.5p . Suppose that v$ is continuous, 
bounded and nonnegative, and fulfils 

- c\\x\~ l <vq(x) <\x\~v- - c 2 \x\~ l for \x\ > 1, (2.6) 

where I is as in (|1.6p and c±,c 2 > 0. Assume also that vq(x) < \x\~^ for all x ^ 0. Let v 
denote the solution of (|2.4p . Then: 

(i) There exist K\ , K 2 > such that for t > 1 we have 

Kx < \\v(;t)\\oo < W, 7 = 7 (i) = ^-M~2)(n-^ (2J) 

^mJ For each ro > one can find C\,C 2 > suc/i that for t > 1 and |x| > ro i/ie following 
holds 

Cie~ xt < \x\-» - v(x,t) <C 2 e~ xt , X = \(l) = (I - fi - 2)(n - I). (2.8) 



Let us comment on the contents and scope of the result. 

1. First of all, it states the two main aspects of the convergence of the solution v(-,t) 
towards the singular steady state Vq\ (|2.8|) establishes the uniform convergence of v(-,t) 
towards Vo in the complement of a ball centered at the origin, with a precise rate that 
depends explicitly on the tail decay exponent /. On the other hand, estimate (|2.T|) gives 
the exact rate of growth of the solutions as t — > oo to account for the approach to the 
singular value Vb(0) = +oo. 

2. An important feature of the result is the existence of a continuum of grow-up rates 
for ||v(-,t)|| co , and a corresponding continuum of stabilization rates of v(-,t) towards Vq 
in the outer region. Note furthermore that as I approaches the lower value jjL + 2, the 
rates go to zero. This limit case is on the other hand easier and does not produce any 
convergence, since we can consider the example of the generalized Barenblatt solutions Vd 
given in (|2.3j) . Indeed, they satisfy < Vd < Vo and 

V (x) - V D (x) = C\x\~^ + V + o (\x\~^ +2 ^ as \x\ -> oo. 

Since they are stationary, no convergence to Vo holds in this case. 

3. The conditions on I imply that the perturbation Vq — vq is never integrable, contrary 
to the usual assumptions made in variational methods. Let us now examine the maximal 
grow-up rate that we have achieved. Note first that j{fi + 2) = j(n) = 0. The maximum 
of 7 in (I2.7p is attained at I = L, and 

j(L) =ii(n + 2-/Jt- 2^/2(n- ^ . 

This is lower than the maximal growth rate of any bounded solution that is given by 
the growth of the spatially homogeneous solution v(t) = ce^ nt . We conjecture that j(L) 
is the largest exponent that can be achieved by the solutions under the conditions of the 
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theorem, even if we allow I to be larger than L. The bound from below follows immediately 
from the lower bound in ()2.7|) . but to obtain the corresponding bound from above is still 
an open problem. 

4. As m — > m* we have /i + 2 — > n and the interval (fi + 2, L] shrinks to the empty set 
while the admissible values of the exponents 7 and A go to zero. 

5. Results similar as in Theorem 12.11 were obtained for the standard Fujita equation 

( N ~ 2 ) 2 - 4iV + 8 v / A r3 T 
u t = Au + U P, x G R , n>10, p>± { N_ 2 )(N-W) ' 

in [TlEli. 

6. Finally, we apply the results of Theorem 12.11 (i) to prove Theorem II. 11 Notice that 
under the assumptions of Theorem 1 if we take the prescribed value of T then vo satisfies 
the hypotheses of Theorem 12. 1\ so that the solution v is global in time and stabilizes to 
Vo; this means that the extinction time of u is precisely T. The extinction rate of u is 
obtained by rewriting the bounds in (12. 7\) . Recall that 



IK,t)||«, =R(Tr n \\v(;t)\\ 00 ~ (T-T)"^, 

and T — t = Te~ 2 ^ n ~^ 1 . The conclusion follows. 

3 Auxiliary results for the rescaled problem 

After the previous transformation, in the radially symmetric case we end up with the 
problem 



Vv :=v t -± ({v m ) rT + ^{v m ) r ) - firv r - pmv = 0, r > 0, t > 0, 
i/(r, 0) = Vo(r), r > 0. 



(3.1) 



An important role is played by the quadratic equation 

a 2 - (n - 2 - fi - k)q + 2k = 0, (3.2) 

where 

K = ('-"-'>("-'> (3.3) 

I — fl 

is positive if [J, + 2 < I < n. The roots a_ and a + of (|3.2p are given by 



n — 2 — a — k± \/(n — 2 — u — k) 2 — 8k 
a± = ^ *\ ^ >- , (3.4) 



and the following way to rewrite a± indicates why the value I = fi + \/2(n — fi) plays an 
important role in the sequel (cf. Section [5]). 
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Lemma 3.1 Assume (|1.5p . and that I £ (li + 2,n). T/ien i/ie roots a± o/ (|3.2p can 6e 

expressed as follows: 

i) If fi + 2 < I < ll+ \j2(n - li) then 

2(n-l) 

a_ = I — li — 2 and a+ = — . 

I — ll 



ii) If n+ \/2(n — fx) < I < n then 



2(n-l) 

a„ = and a+ = I — u — 2. 

I — li 

Proof. Since it can easily be checked that a± := I — li — 2 < 2 ^"~^ =: a 2 if and only if 

I < fj, + \j2{n — li), we only need to check that both a\ and a 2 solve (|3.2p . As to a±, this 
follows from 

a\ - (n - 2 - 11 - k) + 2k = (I - ll - 2) 2 - (n - 2 - ll - k){1 - ll - 2) + 2k 

= (I - ll - 2){l - li - 2 - (n - 2 - LL - K)\ + 2k 
= (I - Li-2)(l -n) + (I - Ll)K 

and the fact that (I — li)k = (I — li — 2){n — I). Using a 2 + 2 = 2 ^_^ , we moreover compute 
a\~ {n-2- li- K)a 2 + 2K = (-^ '-) - {n-2- li)— + (a + 2)k 



2(n-l) 



(2{n - I) - (n - 2 - li)(1 - li) + (n - ll)(1 - li - 2)) , 



{l-Lif 

from which we immediately find that also a 2 solves f)3. 2 [) . □ 

The following two lemmata apply to parameters n, m and k more general than required 



in §\J$ and (J3J3J). 

Lemma 3.2 Let n > 1, m > 0, k > and o~q > 0, and set 

a(t) := a e^ Kt , t>0, 

and 

£{r,t) := a»(t)r, r > 0, t > 0. 

Suppose that tp : [0, 00) — >■ [0,oo) is twice continuously differentiable in (£o,£i) with some 
£0 one? £1 satisfying < £0 < £i- T/ien /or 

u(r, t) := a(t) (> 2 (r, t) + r/>(£(r, t))) 2 , r > 0, t > 0, 
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we have the identity 

Vv(r, t) = (e 2 (r, t) + ^(£(r, t))) " # ^((£(r, t)) (3.5) 

/or a// (r,t) G 5 := {(p,r) G (0,oo) 2 | £0o,r) G (Co, 6)}, 

^(0 := (e 2 + ^) + + 2^ - (/x + - (3.6) 

/ore G (Co, 6)- 

i --1 1 

Proof. Using £ t = -av otr = and at = una, we compute 

v, = n(e + m)~ f ' -f<'(« 2 + *(£)))~ S ~ 1 (2£ + ^(«))6 

. ff ff ( { . + ^ t)) )-t-'{^(«. + ^)_Is (a{+W0 ) { } 

_ i£ _ 1 

= |ff(e 2 +^(0)) a {2«^(o-«e^(o} (3-7) 

whenever (r, i) G 5. Since £ r = at 1 , we moreover have 
as well as 

(v"V = + (2 + fe(fl) 

+=(= + iy- + i (e + «f)))^"(2{ + 

at such points. Thus, in view of the identity m+ |- = 1 and, equivalently, ^ = ^ — 1, we 
find that 

^{( Vm )rr + ^(v m )r) = ~\ ^ + ^)) ^ (2 + %(£)) 

+(f) 2£j (^ 2 +^)) ~ f_1 ( 2 ^+^(o) 2 

+| + ^(o) 2 - ^ (e + ^(0) (2e + mo) } 
= f^+^e))"^ 1 !- (e 2 +^(o)(^(o + ^^(o) 

_ 2 ( n _ ^)e 2 - 2nV(C) + 2^(0 + f (£)} (3-8) 
8 



if (r, t) G 5. As r£ r = £, we finally have 

firv r = /U r|-^(e 2 + v(o) _f_1 (2e + ^(o)er} 

and therefore obtain from (13.71) and (13.81) that 



If 

+(e 2 + m) + ^^(o) + % - # 

+2nV(C) - 2^(0 - f ^(0 + 2^ 2 + /i^(e) - 2< 2 - 2n^(e)}, 
which after a straightforward rearrangement yields (|3.5p . □ 
Lemma 3.3 Let n > 1, m > and e > 0. XTien 

^):=i-ee 2 , e>o, 

^(£) = 2(k - ne) - 2(n - /i)e(l - e)£ 2 /or a« £ > 0, (3.9) 
where A is defined by (|3.6p . 



Proof. We directly compute 

AiP(0 = (£ 2 + 1 + ^ 2 )(-2e - (n - l)2e) + 2«(1 - ef) + (/x + k)2^ 2 - |(-2e£) 2 
(1 - e)£ 2 + l) (-2ne) + 2k - 2«e£ 2 + 2k^ 2 + 2/xe^ 2 - 2/ue 2 £ 2 
= -2ne(l - e)f - 2ne + 2k + 2/ze£ 2 - 2 / ue 2 ^ 2 for £ > 0, 
and thereby immediately obtain (|3.9|) . □ 



4 Lower bound 

Once we have the preparatory material, we proceed next to establish the lower bound for 
the solutions mentioned in Theorem 12. 1[ This is the content of Proposition 14.21 Section [5] 
will contain the proof of the corresponding upper bound, Proposition! 



Lemma 4.1 Suppose that condition (jl.5p on n and m holds, and that Z G (zi + 2,n) and 
a G [a_,a + ] with a± given by (13. 4p and k as in (13.3j) . Then, there exist a > 0, £o > 
and a positive tp G C°([0,oo)) n C 2 ([0,oo) \ {£o}) satisfying 

Ail><Q in (0,oo)\{£ } (4.1) 
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and 



lim inf ?/>(£) > lim sup ip (£) , (4.2) 



and such that 

-> a as £ ^ oo. (4.3) 
i?ere ; .A is i/te operator defined in (|3.6p . 



Proof. We let e := - and fix a > small such that 



— +e — >a«+2 < 1. (4.4) 



a/ \2e/ 
Then the function ip defined by 

has a unique local minimum at the point at which = -aa^~ a_1 + 2e£ = 0, that is, 

at the point 

aa \ ST2 



S,min 



with corresponding minimum value 



- .(=)"* - 1 + <(f )* - {(I)* + 

by (|4.4p . Therefore, 

£ := inf{£ > | v?(0 < 0} 



OQ +2 - 1 < 



lies in (0,£ m i n ) and we have 
Accordingly, 



Vtfa) < 0- ( 4 - 5 ) 



if^G[0,&], 
if £ > Co, 



defines a positive continuous function on [0,oo) which, by (|4.5p . satisfies (|4.2p and clearly 
also fulfils (|4.3p . Moreover, using / > /i + 2 and (|3.3p we find 



(l-/j,)(n-K) = (I - fi)n- (I - fi-2)(n - I) = I 2 - (fi + 2)l + 2n 

+ 
2 



2ra > 0, 
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hence s = ^ < 1. Thus, recalling Lemma 13, 31 and the fact that /i < n, we obtain 

Aip(0 = -2(n - n)e(l - e)f < for all £ G (0, Co)- (4.6) 

As to large £, we compute 

Ajj)(£) = (i 2 + aC a ) (a(a + l)aC a ~ 2 - (n - l)oa£ _ 

/' 
— i 

2 



a(a + 1) — (n — l)a + 2k — (/i + k)o ja£ a 
+{a(a + 1) - (n - l)a - |a 2 }a 2 r 2a " 2 

a 2 — (n — 2 — // — «;)a + 2fi;|a£~ Q 

.|£Z^ a 2 + („_2)a}a 2 r 2a " 2 for^>eo- 



Since a 2 — (n — 2 — fi — «)a + 2k < due to our assumption a £ [a_,a+], and since 
\i = > 2 and n > 2, from this we immediately infer that 

A±{£) < for all f > Co, 

which combined with (|4.6p proves (|4.ip . □ 

Proposition 4.2 Assume again condition (jl.5p . Suppose that vq G C°([0, oo)) is positive 
and fulfils 

Vo(r) > r _/x — cor - ' /or all r > 1 (4-7) 

with some cq > and Z as in (jl.6p , XeZ i> denote the solution of (|3.ip . Then: 
(i) There exists C\ > suc/i iZiat 

fi(l — f.L — 2)(n — l) , 

v(Q, t) > Ci e «-#« /or aZZ t > 0. (4.8) 
(mJ For eacZi ro > one can /md C2 > fulfilling 

v(r, t) > r~" - C 2 e -('-**~2)(n-Qt y or a // r > rQ am f t > . (4.9) 

Proof. (i) As / is as in (|1.6|) . we can apply Lemma 14. II to a_ = 1 — ^,-2 and obtain 
£0 > and a function ^ with the properties provided by that lemma. Since vq is positive, 
there exists c\ > such that 

v (r) > ci for all r G [0,1]. (4.10) 
Let us pick C2 > large such that with Co as in (|4.7p we have 

(l + c 2 )"t < l-c , (4.11) 
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and then choose £ > small fulfilling 

rV(0>C2 foralUe(0,O, (4.12) 
which is possible since ip is positive on [0, oo). We next define 

z := max C 2 ^(0 (4-13) 
fe[|,oo) - 

and let C3 > 0, C4 > and C5 > be small enough such that 

(1 + z)~$ <l-c 3 z for all z £ [0, z ] (4.14) 

and 

£ 2 + ±(0 > c 4 for all £ > (4.15) 

as well as 

</>(£) > c 5 f" for all £ > £ (4.16) 

where we make use of (|4.3p and, again, the positivity of ip. Finally, we take a small number 
do > satisfying 

<x < cicf (4.17) 

and 

aQ <(^)^ (4.18) 
V c / 

and define 

a(t) := aoef^, t > 0, 
with k = H^Ez^HzH as i n (|3.3j) . Then 

«(r, t) := a (t) (V(r, t) + ^(£(r, t))) 2 , r > 0, i > 0, 

with £(r, i) := a»(t)r, is continuous in [0, oo) 2 and smooth at each point (r,t) G [0, oo) 2 

_ n_ 

where r 7^ ?"o(i) := ^o" 7 ''(O- An application of Lemma 13.21 and Lemma 14.11 shows that 
with A as in (|3.6I) , 

Vv = ~<r(i) 2 (r, i) + ^>(£(r, i))J 2 -4V>(£(r, i)) < whenever r / r (t), 



which implies that v is a subsolution of (|3.ip in the Nagumo sense, because from (|4.2p we 
infer that 

limsup?; r (r, t) < liminf v r (r, t) for all t > 0. 

r/V (t) r\r (t) 

Accordingly, if we can show that u does not exceed v initially then the comparison principle 
will tell us that 

v>v in[0,oo) 2 (4.19) 
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and, in particular, 

v(0, t) > v(0, t) = a{t) = a e^ Kt for all t > 0, 
which will yield (|4.8p in view of the definition of k. It thus remains to show that 

v(r,0) < v (r) for all r > 0. (4.20) 
To this end, we first consider the case r < 1, when (I4.15p . (I4.17P and (|4.10p imply 

v(r,0) = <Jo(f(r,0)+±((i(r,0))) * 

< (Tqc^ < ci < Vo(r) for all r G [0, 1]. (4.21) 

Next, if r > 1 is such that r > |cj " then £(r, 0) = a M r > |, and hence from (|4TT5|) . (gH]) 
and (|4.16p we obtain 

v(r,0) = r-^l + C\r,0)l(d(r,0))) 2 < r~" (l - c 3 C 2 (r, 0)Y>(£(r, 0))) 
< r~" (l - c 3 c 5 r Q " 2 (r-, 0)) = r"" - cs^^r"^" 2 , 



which in view of our choice a = I — fx — 2, (|4.18p and (|4.7p gives 



l — fi 



v(r, 0) < r M - c 3 C5(T " r ' < r M - c r 1 

< v (r) for all r > minjl, |ctq 77 1 . (4.22) 

Finally, if r > 1 is such that r < £c M then we use (|4.12|) and (14. lip to see that 

v(r,0) = r^(l + r 2 M)^(£(r,0))) 2 < r^(l + c 2 )~t 

_ j_ 

< r -A *(l - c ) for all r < |cx " , 

whereas by (|4.7p . 

«o(r) > r~ M (l - c r~ ( ^ )N ) > r"^(l - c ) for all r > 1, 



so that v(r,0) < v (r) also holds if 1 < r < £ct M . Together with P~2T|) and (|4~22|) this 
establishes (^20|) . 

(ii) To see (|4.9p . we observe that in view of (|4,3p there exists c@ > such that 

V>(6 < c 6 r° for all £ > 0, 
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where still a = a_ = I — \i — 2. Then by (|4.19p and the convexity of z \-± (1 + z) 2 for 
z>0we obtain 



v(r,t) > v(r,t)=r-'*(l + S- 2 ±(Sj) ' > ^ ~ |r _ TV(0 

= r"^ - ^c 6 a~^r- 1 e -(J-**-2)(n-9* for all r > and t > 0. 
Given ro > 0, this easily yields (|4.9p upon an obvious choice of C%. □ 

5 Upper bound and proof of Theorem 12.11 

Lemma 5.1 Assume (jl.5p . and let I € (/i + 2,n) and a_ fre as defined in (j3.4|) irot/i /c 
given fry (|3.3|) . T/ien i/iere exisi /3 > a_ and > wi/i the following property: Suppose 
that A > and > are such that 

B a-+2 

and let 

f3B 

Cl == V 

T/ien t/ie function tp^t defined by 



<(wR) ~ - Bf" S , f > 0. (5-3) 

satisfies 

WwM£l)=0 and &ouM0<0 for all O^, (5.4) 
and moreover we have 

Moutit) > for all £ > 6 (5.5) 

with A given by (|3.6p . 



Proof. Throughout the proof, let us abbreviate a := a_ for convenience. Then since 
I G (/t + 2, n), we have < a < a + , and 

p(/3) := /3 2 - (n - 2 - /t - k)(3 + 2k, £ R, 

has the properties 

p(a) = and p(/9) < for all /3 G (a, a+). (5.6) 

As furthermore 

q(j3) := 13(13 + 2 - n) - a (a + 2 - n) + //a/3, /3 G R, 
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evidently satisfies q(a) = fia 2 > 0, by continuity we can choose some /3 > a such that 

< 2a + 2, p(fi) < and q(fi) > 0. (5.7) 
With this value of /3 fixed, we pick Cp > large such that 



and 

n /3+2 i 1 

2C3 /a\fr^/J_ 4 

where 



< 1 (5.8) 



ci := -p(/3), c 2 := ^y^a 2 + (n - 2)a and c 3 := ^y-V + (n - 2)/9 (5.10) 

are all positive according to (|5.7p and the inequalities = > ^ anc ^ fi > 2. Then, 
given A > and > fulfilling (|5.ip . we let ip out be defined by f)5.3j) and compute 

(VwMO = -a^r^'+^r^ 1 and 

(Vw)tt(f) = aia + lMr-'-W + ljBr 3 " 2 (5.H) 

for £ > 0. From this it can easily be deduced that in fact i> out attains its maximum at 
£ = £i and decreases on (£i,oo), where £i is as in (|5,2p . Using (|5.1ip we furthermore 
obtain 



MoutiO = {e + AC a -BC l3 ){a{a + l)AC a ~ 2 -^ + l)BC p - 2 

-{n - l)aAC a ~ 2 + (n - l)(3BC P ~ 2 } 
+2kAC° ~ 2kBC P - (p + n)aAC a + (p + k)I3BC P 
_| (a 2 A 2 r a - 2 ~ 2af3ABC a ~' 3 - 2 + P 2 B 2 C 2 ^ 2 ) 

= {a(a + l) - (n-l)a + 2K- (p + K)a}A£~ a 

+{ - P(p + 1) + (n - 1)13) - 2k + (/x + k)^}^-' 9 
+| a ( Q + 1) - ( n - l) a - |a 2 }A 2 r 2a ~ 2 
+{/3(/3 + l)-(n-l)/3-^ 2 }i? 2 r 2/3 - 2 

+{ - a(a + 1) + (n - l)a + /3(/3 + 1) - (n - l)/3 + /ia^ABr -^ 2 

= p(a)^r a -P(P)BC P ~ {^« 2 + (" " 2)a}^ 2 r 2a ~ 2 

-{^yV + (n - 2)/3}s 2 r 2/3-2 + g^ABr """ 2 for all £ > 0. 
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Now (|5.6p and (|5.7p imply that the first term on the right vanishes and that the last 
is nonnegative, because A and B are positive. Hence, recalling (I5.10P we arrive at the 
inequality 

A^ out (0 > Cl BC P ~ c 2 A 2 r 2a ~ 2 - c 3 B 2 r W ~ 2 for all £ > 0. (5.12) 
Now if £ > £o then (|5.8p along with (|5.ip and our restriction (3 < 2a + 2 ensures that 



Moreover, for such £ we find 

c 3 B 2 tj~ 2 ^ 2 2c3 £-2 < 2£3 B ((W\ ifer 
±ciB^P ci - ci VaA/ 



by (HlH]). Together with (15TT3J) and (157L21) . this shows that indeed ^ out > for all £ > £i, 
as claimed. □ 

Lemma 5.2 Suppose that (jl.5p holds. Let I G (/i + 2,n) and a_ &e as m (|3.4p urai/i At 
gfiwen 6?/ (|3.3p . T/ien i/iere exist ^4 > 0, £i > and a positive function i/j G C°([0,oo)) n 
C 2 ([0, oo) \ {£i}) such that with A as in (13.61) . 

A0>O in (0,oo) \ {£1} (5.14) 

and 

lim sup t/jt (£) < lim inf tp* (£) (5.15) 

as u>e// as 

f -^(0 -»• A as £ -)■ oo. (5.16) 

Proof. Again we write a := a_ for simplicity. Since jjL < n, it is possible to fix 
ci G (0, 1) such that 

c ' £ (5 - 17) 

and since I G (/i + 2,n), there exist /3 > a and (7g > such that the conclusion of 
Lemma |5 . 1 1 holds . We now define 
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which is positive because < a < (3, and 



Next, we pick A > large fulfilling 



^>{^c?Q^KTO} f (5.20) 



and 



and let 



and 



as well as 



so that 



A> [CpK^r)^^ (5.21) 



B :=(-)" (5.22) 



(5.23) 



\aAJ 



(5.24) 
(5.25) 



Then the function ip : [0, oo) —> R given by 

1 Vw(£) = AC a - bc p if ee (£i,oo), 

is continuous on [0, oo), because (|5.25p ensures that 

^in(Ci) = 1 ~4, 

whereas invoking (|5.18|) . (15.22j) and (|5. 19j) we find 

Vw(6) = - S(g)" 7 ^ = c 2 ^"^ = c 2 i^ = l-d 

Note that since ci < 1, this also implies that ip is positive on [0, oo). Next, from Lemma f3.3l 
and (|5.25p we obtain 

A$ in (0 = 2( K -ne)-2(n-fi)e(l-e)e 
> 2(k - ne) - 2(n - n)e^\ 
= 2(k - ne) - 2(n - n)<% for all £ e (0, £i), 
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where by (|5T23l) . (|5T22l) and flOUD . 



e = c\\ 1 =cf(^)"""K^^),4"^ < 



2 

/I \ P-a , , 2 



Hence, 

>2(/s-|)-2(n-/x)c?>0 for all (0,6) (5.26) 
according to (j5. 17|) . 

Now the requirement (|5,2ip along with (|5.22j) guarantees that 

B a+2 ^_«+2 2 (ff-a) 

so that Lemma 15. II becomes applicable to tell us that 

A^ out (0 > for all f > £i (5.27) 

as well as 

= 0- (5-28) 

As a consequence of (|5.26p and (|5.27p . we see that (|5.14p holds, while (|5.28p combined 
with the fact that 

(^in)e(Ci) = -2e£i 2 < 
yields (|5.15p . The assertion (|5.16p immediately results from the definition of ijj. □ 



Proposition 5.3 Suppose that (jl.5p holds, and that v is the solution of (|3.ip . where the 
initial data vq S C°([0, oo)) are nonnegative and such that there exist I as in (jl.6p and 
c\ > fulfilling 

v {r) < - cir~ z /or r > 1, (5.29) 
and which in addition satisfies 

v (r) < /or aZZ r > 0. (5.30) 

There exists C\ > suc/i i/ioi 

fj,(l — fi~~2)(n — l) , 

v(r,t)<Cie for allr>0 andt>0. (5.31) 

iij For a// ro > i/iere exists C2 > u>ii/i £/ie property 

v(r, t) < r~v - C 2 e -( i -^- 2 )( n - / )* / or tf r > rQ andi > 0. (5.32) 
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Proof. i) Since I is as in (jl.6p . the number a := a_ satisfies a = I — fi — 2 by 
Lemma 13.11 Hence, applying Lemma 15.21 we find £i > and a positive function tp £ 
C°([0,oo)) n C 2 ([0,oo) \ {£].}) with the properties (15TLH) and ([57L5J) and such that 

< c 2 for all £ > (5.33) 

with some c 2 > 0. Taking c 3 > large such that 

vo(r) < c 3 for all r > 0, (5.34) 

we can find tq > small enough fulfilling 

r < (2c 3 )~» (5.35) 
and then, by (|5.29p and (|5.30p . fix C4 6 (0, c\] such that 

v o( r ) < r_/i ~~ C4r~' for all r > ro. (5.36) 
We pick £ > and C5 > sufficiently large satisfying 

£>(//c 2 )^ (5.37) 

and 

e 2 + ^(0<cg for all CG [0,|] (5.38) 
and finally choose a large number ctq > with 



and 



We now define 



(Tq > c 3 c 5 2 (5.39) 



..>(f|)- (5.40) 



o-(t) := cjoe^*, t > 0, 

and 

v(r, t) := a(t) ^(r, t) + ^(£(r, t))) 2 , r > 0, t > 0, 

again with £(r, i) := a^{t)r, and claim that 

v(r,0) > v (r) forallr>0. (5.41) 

Indeed, if r < |<r " then £(r, 0) < | and hence (f5738]) . ([09|) and ([5341) imply that 

"(r,0) = a (^ 2 (r,0) + VK£(r, 0))) ~ 2 > a c~ * > c 3 > «a(r) for r < £a~*. (5.42) 
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Next, in the case when £ct m < r < ro we have £(r, 0) > £, so that using the convexity of 
0<m(l| z)~2 along with (15T331) . (IOTP . (1Q51) and (|534"]l . we can estimate 

tf(r,0) = r-"(l+r 2 ('',0)?(«r,0))) 1 > r""(l - |r 2 (>-. <»«(«>•, 0))) 

> ^ > c 3 > «o(r) if la" 77 < r < r . (5.43) 

Finally, for r > £<7 m fulfilling r > ro, by the same convexity argument in conjunction 
with the fact that a = I — /i — 2, from (|5,40p and (|5.36p we have 

v(r,0) > r-^(l-^c 2 r a - 2 (r,0))=r-^-^c 2 a^r- 1 



c^r 1 for r > max jYo, £<7 " \. 



Together with (|5.42p and (|5.43p . this proves (|5.4ip . Since by Lemma [3721 and Lemma l3~3j 
recalling (|3.6p we have 



Vv = ^a(t)(^ 2 (r,t) +i>(£(r,t)f) 2 Atjj(r,t)>0 whenever r ^ n(t) := £i<r~»(t) 
and 

liminf v r (r, t) > lira, sup v r (r,t) for all t > 
r/n(t) r \n(t) 

according to (|5,15p . the comparison principle applies to yield 

v > v for all r > and t > 0. (5.44) 

This immediately leads to (|5.3ip . 

ii) To obtain (|5.32p . we fix ro > and first pick C5 > small enough fulfilling 

(1 + z)~2 < 1 - c 5 z far all are [0,1], (5.45) 
and then fix to > large such that 



I := c 2 a " r~ ( '~ M) e -(i-A«-2)(n-i)to < 1? ( 5 . 46 ) 



where c 2 and 00 are as determined by (|5.33p . (I5.39P and (|5.40p . Then for all r > tq and 
t >to, still writing a = a_ = I — fi — 2 we have 

C2 r Q ~ 2 M) < c 2 C a ~ 2 (r ,t ) = I < 1 
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and thus 



by (|5.45p . Therefore, (|5.44p and (|5.33[) entail that for such r and t we have 

v(r,t) < v(r,t) =r'^(l + C 2 ^(0) ' 

< r-»{l + c 2 r a ~ 2 )"^ < r-"(l - c 2 c 5 r a " 2 ) 

= r"^ - c 2 C5CT"^r-^- a - 2 = r - " - c 2 c 5 a~^ r' 1 e -('-f-2)(«-0*. 

This shows that (|5.32p is valid for some sufficiently large C 2 > 0. □ 
As we have said, Propositions 14,21 and 15.31 together imply Theorem 12.11 

6 Comments and open problems 

1. The construction of the new extinction rates for m £ [m*, m c ) is open. The relevance of 
m* in the asymptotic analysis of stability of the Barenblatt solutions has been documented 
in EES H| . 

2. We have not performed the analysis of positive perturbations of the tail of the singular 
solution Vo- Preliminary calculations show that we can have in that case global grow-up 
if the perturbation is large, i.e., if I — fi > is small. The case I = jjL is explicit; indeed, it 
is easy to check that the solution with initial value vq(x) = A \x\~^ is 

v(x,t) = (Ce 2{n -^ + l) 1/(1 - m) \x\-^, C = A 1 ~ m -l. (6.1) 

For A > 1 this solution blows up everywhere as t — > oo with rate 0(e 2 ^ n ~^ t ), while for 
A < 1 it vanishes in finite time. 

3. The analysis of perturbations of the Barenblatt profiles, Vd with D > 0, with large tails 
of the form vq(x) — Vd{x) = 0(\x\~ l ), is an interesting related problem. The difference with 
the above analysis is that the ^-profile is regular, so no grow-up is expected if I > fj, + 2. 
Since the behaviour of Vd at infinity is similar to the singular one, Vo, and Vo is still 
stationary, we also expect a continuum of convergence rates depending on I from a certain 
range. In this case we have to mention that for I > n there is a variational theory developed 
in the recent papers [21 El 0] that proves convergence with rate using the techniques of 
entropies, linearization and functional inequalities. 

4. We could have used another of the possible scaling options, which is not adapted to 
the Barenblatt profiles but is still adapted to the singular solution. The simplest choice is 

w(y, s) = [(1 - m)(T - T )]-V(i— ) > n(y , r ), s = (1 - m) log[(T - r)/T] = (2/(3) t, (6.2) 
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which leads to the equation 

^ = A(w m /m)+w. (6.3) 

Putting w m = Z and p = 1/m we get a variation of the Fujita equation 

dZ p 

^— = aAZ + bZ p . (6.4) 

OS 

Studying this equation is equivalent to the study of the v equation. It is interesting to 
translate the results we have obtained and to compare with the standard Fujita equation 
ut = Au + u p . 

5. Our methods are not variational and our solutions do not belong to the usual spaces 
of that theory, like spaces of finite relative energy or finite relative mass. 
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